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Abstract 

In this paper we obtain explicit expressions for tau-functions related to Picard 
type solutions of the Painleve VI equation in terms of theta functions and their 
derivatives. 



1 Introduction 

In this paper we study a special case of the Painleve VI equation [1,2] 

^"(t) = l(4r.+ -71^ + -71^) ^'(ty -(7 + tAt + ) <l'it) + 



2vg(t) qit)-i q{t)-tj \t (t-l) q{t)-t 



q{t){q{t) - l){q{t) - t) 

when 



^ t t-l , t{t-l) 
+ P^T7;+lT77 77; + 5-, 



a = 0, /3 = 0, 7 = 0, 5=^. (2) 

This case was originally considered by Picard [3]. Due to a special choice of parameters 
(2) a general solution of (1) is known 

go(i) = p{,CiUJi + C2CU2; tui, ^^2) + (3) 
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where p{u; oui, 002) is the Weierstrass eUiptic function with half-periods 001^2, Ci,2 are com- 
plex constants and uji^2{t) are two linearly independent solutions of the hypergeometric 
equations 

t{l - t)uj"{t) + (1 - 2t)uj\t) - ^uj{t) = 0. (4) 

The properties of the Picard solutions have been studied recently by M. Mazzocco [4]. 
In particular, she investigated its monodromy properties and algebraic solutions which 
correspond to Ci and C2 being rational numbers. 

Algebraic solutions of Painleve VI play an important role in many applications in 
theoretical physics (see, for example, [6,7]). In such cases a calculation of related tau- 
functions can be simpler due to a presence of an algebraic relation between the solution 
q{t) and the variable t. 

The goal of this paper is different. We aim to present explicit expressions for tau- 
functions related to the Picard solutions and its images under birational canonical trans- 
formations [5] for generic values of complex parameters ci and C2. To our knowledge this 
has not been done before. 

More explicitly, we consider a sequence of tau-functions [5] obtained by a parallel shift 
/a from the Picard solution. We calculate the first two tau-functions and other members 
of the sequence can be obtained using the standard Toda-type second order relations. 



2 Properties of the Painleve VI equation 

In this section we briefly review the main properties of the equation (1) which we denote 
as Py/(Q;,/3,7,(5). 

Following [5] one can introduce two different parameterizations of parameters in the 
Painleve VI equation: {hqiH\i k^o, 0) 

and {hiMMM) 

Ko^bi + b2, Ki = 61 - 62, Koo^h-h, 6 ^b3 + b4 + l. (6) 

This equation is equivalent to the Hamiltonian system Hvi{t;q,p) described by the 
equations 

dq dH dp dH . , 

dt^~d^' dt^ ~~d^ ' ^ ^ 

with the Hamiltonian function 

Hvi{t; q,p) = i)p'^- 
-{ko {q -l){q-t) + K, q{q - t) + {9 - l)q{q - l)}p +K{q-t)], (8) 
where q = q{t), p = p{t) and 

K^^{Ko + Ki + e-lf-^Kl. (9) 
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One can introduce an auxiliary Hamiltonian h{t), 

h{t) = t{t - i)H{t) + 62(61, 63, h) t - ^62(61, 62, h, 64), (10) 

where ei{xi, . . . , Xn) is the i-th elementary symmetric function in n variables and a set of 
Xj's can be a subset of biS as in (10). 

Okamoto [5] showed that for each pair {q{t),p{t)} satisfying (7), the function h{t) 
solves the Ey/ equation which is 
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h'(t)^t(l-t)h"{t)j +^h'(t)[2h(t)-{2t-l)h'(t)]+bib2bsb^j =n(^'(^) + ^^) (^^) 

k=l 

and q{t) solves Pi//(a, /9, 7, 5). 

Conversely, for each solution h(t) of (11), such that -^hij:) 7^ 0, there exists a solu- 
tion {q{t),p{t)} of (7), where q{t) solves (1). An explicit correspondence between three 
sets {q{t), q'{t)}, {q{t),p{t)} and {h{t), h'{t), h"{t)} is given by birational transformations, 
which can be found in [5]. 

The group of Backlund transformations of Py/ is isomorphic to the affine Weyl group 
of the type F4: Wa{FA). It contains the following transformations of parameters (only five 
of them are independent) 

: 61 -H- 62, ^2 : 62 -H- 63, ^3 : 63 -H- 64, W4 : 63 ^ -63, 64 -64, (12) 

: Kq Ki, x^ : Ko k^, x^ : Kq 9 (13) 
and the parallel transformation 

h:h= (61, 62, 63, h) ^ b+ = (61, 62, bs + 1, 64). (14) 
The auxiliary function h^{t) corresponding to parameters b"*" — /3(b) is given in [5] 

h+{t) = h{t) - q{q - l)p + (61 + 64)? - ^{bi + 62 + 64). (15) 

Following [5] one can calculate h+{t) in terms of h{t) and its first and second derivatives 

, , , _ t{t - l)h"{t) + 2h{t)[b3{bs + 1) + h\t)] + 63(1 - 2t)h'{t) - 616264 
' " 2{h'{t) + 61) 

and vice versa 

^^^^ _ t{t - l)h':,{t) + 2h+{t)[b,{b, + 1) + h'4t)] - (63 + !)(!- 2t)h'4t) + 616264 

2(/.V(t) + (63 + 1)2) ■ ^''^ 

For each solution of the Pyi equation one can introduce a corresponding tau-function via 

H{t,q{t), p{ty,h) = j^logT{t,h). (18) 



Obviously tau-functions are defined up to an arbitrary normalization factor. 
Following [5] let us introduce a family of tau-functions T^(i) 



T„(i) = exp<^ / dtH{t,q{t),p{ty,K 



where 



l^{h) = {bi,b2,h + m,bi), meZ. 



(19) 



(20) 



As shown in [5] they satisfy the second order Toda-type equation 



d 
di 



t{t- l)j^\ogTm{t) 



+ (b, + bs + m){bs + b, + m) = c{m) T"-+i(|^)Tjn-i(t) ^31) 



where c(m) is a nonzero constant. 



3 Elliptic functions and useful identities 

In this section we list all definitions and properties of elliptic functions used later in the 
text. Following [8] we will use the standard theta functions 6i{x\T) with quasi-periods tt 
and TTT. 

The elliptic modulus k and its complement k' are defined in a standard way by 

It will be more convenient to use the parameter t = k"^ as the second argument of elliptic 
functions and hereafter we will follow this notation (except for the theta- functions), i.e. 

r = z|^, ? = e- m^^elm, K'{t)^K{l-t) (23) 

where K{t) and K'{t) are the complete eUiptic integrals of the first kind of the parameters 
t and 1 — t. 

We introduce Jacobi elliptic functions 

where 

„=H^« (25) 

TT 

and define the fundamental elliptic integral of the second kind [8] by 

£{u,t) = / dn^ix, t)dx. (26) 
Jo 

It satisfies 

S{K{t),t) = E{t), (27) 
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where E{t) is the complete eUiptic integral of the second kind. 

Using (23) and identities for complete elliptic integrals of the first and second kind 

W) = ^-M. 9,SW = M^ (28) 



2t{l-t) 2t ' ' ' 2t 

E{t)K{l -t) + E{1 - t)K{t) - K{1 - t)K{t) = 7r/2, 



one can obtain 



dtrit) 



ITT 



(29) 



(30) 



4t{t - l)K^{t)' 

We also need the derivatives of Jacobi elliptic functions with respect to the parameter 
t. Differentiating the formula 

sn(u,t) 

"= / ;7Tr 



dx 



^y{l-X^){l-tX^) 

with respect to t and calculating the remaining integral we get 



(31) 



d , , 
— sn(M,t) 



sii((/. f)cir ((/,. i) cii((/. f)(lii((/,, t) 



2{t - 1) 



2t{t - 1) 



Prom (32) it is easy to obtain 



d 
di 



cn{u, t) 



sn^(-u, t)cn{u, t) sn{u, t)dn{u, t) 



2{t 



2t{t - 1) 



[u(t-l)+S(u,t)] 



and 



d , , , sn'^(u,t)dn(u,t) sn(u,t)cn(u,t) r , ^, ^, 



(32) 



(33) 



(34) 



l(t - 1) 2(t - 1) 

Integrating the well known formula between a logarithmic derivative of 94{x\t) and 
£{u, t) we obtain 



2xK{t)/iT 



2x^ f 
elixir) = e^iOlr) exp ^ -—E{t)K{t) + / E{y, t)dy 



> . 



(35) 



We can use (35) for calculation of the derivatives 9i{x\T)'^ and 94{x\t)'!^. 

Finally, combining (35) with the differential equation satisfied by theta-functions 

Ad 

-^i(M|T) + — ^i(M|T) = 0, i = l,2,3,4, 



m dr 



(36) 



one can calculate the derivative drOi{x\T) in terms of theta-functions and S{u,t). 

We notice that formula (35) is convenient for expansion of 94{x\t) in a series in x up 
to any required order. Say, 



elixir) = ^4(0|r)exp 



2x' 



mm 



(37) 
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4 Picard solution and tau-functions 

The restriction on parameters (2) for the Picard solution (3) can be rewritten in terms of 
parameters 6j (6) as 

6i = 62 = 0, 63 = &4 = -1/2. (38) 

It is convenient to fix a particular branch of the Picard solution (3) by choosing two 
linearly independent solutions of equation (4) 

uJiit) = I 2Fi(i i; 1; = K{t), uj,{t) = i| 2Fi(^, iK'{t) = iK{l - t), 

(39) 

where K{t) and K'{t) are the complete elliptic integrals of the first kind as defined in the 
previous section. 

With a choice of half-periods (39) for the Weierstrass function the expressions for the 
invariants ei, 62, 63 take the following form 

ei = l 3-, 62 = t g-, 63 = g-. (40) 

After some simple calculations we can rewrite the Picard solution of Pyi as 

" sn2(cii^(i) + ic2i^'(i),i)' ^^^^ 

where ci and C2 are the same parameters as in (3). 
Let us make a change of variables 

ci = - + l, C2 = ^ + l, (42) 

TT TT 

where x, y are new parameters. A reason for this is that the resulting tau-functions look 
simpler in x and y. 

With a substitution (42) formula (41) takes the form 



where we defined a new variable z 

2K{t) 



TT 



(x + ry). (44) 



The hamiltonian iJo(^) for the choice of parameters (38) can be calculated from (7-8) in 
terms of go(^) and its first derivative 

^ + (1 - 2t)q,{t) t(t-l)q',(tf 

Atit-l)iqoit)-t)^ Aqoit)iqoit)-l){qoit)-ty ^ ^ 

Using formulas from section 3 one can explicitly calculate the derivative qgit) 
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d 

di 



+ 



sii{z, t)cn{z, t) 



dn\z,t) dn\z,t) 



\oge2{x + Ty\T)]',+ 



2K{t) 



Kit) 



(46) 



Substituting (46) into (45) we produce the following expression for the function Hoit): 



cn^{z, t) 



4{t-l) 4t{t- l)sn'^ {z,t) ' 4t{t-l)' 



+ 



where 



TT 



2K{t) 



[log^i(x + Ty|T)]^ + 



W) 



(47) 



(48) 



and z is defined by (44). 

To calculate the tau-function To(t) for the Picard solution (43) we have to calculate 
the indefinite integral of (47) with respect to the variable t which looks like a hopeless 
problem. 

Now we formulate the central result of this paper. 



Theorem 4.1 The tau-function for the Picard solution (43) is given by 

To{t) = exp{ j Ho{t)dt} = Co{x, y) q^'l^h-^l' 

where CQ{x,y) is an integration constant. 
Proof: 

First we rewrite To{t) as 



(49) 



To{t) = Co{x,y)qy"/-"sn{z,t) 



O^ix + TylT) 



(50) 



The proof is straightforward and reduces to differentiations. We shall do this in a few 
steps. Taking logarithmic derivative of (50) we obtain 



d , ^ / X ^ / X cn(z,t)dii(z,t) 

- log To i = i^dtrit) + ^ ' ; ' 

at TT sn[z,t) 

dt sn{u, t) 



TT 



sn(M, t) 





r) 




r) 



u=x+Ty 



TT 



+ 



dtT{t). (51) 



Using (28-30) one can evaluate the first two terms in (51). The derivative dtsn{u,t) 
was calculated in (32). Differentiating (35) twice and using the equation (36) one can 
evaluate all derivatives in the last term of (51). Combining all contributions and using a 
simple formula 



\og9i{x + Ty\T)l^ log e4{x + ry\r% + 



2K{t) cn{z,t)dn{z,t) 



TT sn(z,t) 

we obtain after simplifications the expression (47) for Ho{t). 



(52) 
I 
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It is quite remarkable that the indefinite integral of the function (47) gives such a 
simple answer (49). We were able to produce this expression by expanding (47) in a series 
in a; at y = and integrating term by term. Analyzing the resulting series we compared 
it to the expansion of 9i{x\t) in x which is similar to the expansion (37). This allowed us 
to arrive at the final answer (49). 

To solve the equation (21) for T^(i) we need to calculate the second tau-function 
Ti{t) corresponding to the solution with parameters 

bi = (0,0,1/2,-1/2). (53) 

First we have from (10) 

ho{t) = t{t - l)Ho{t) + t/A - 1/8 (54) 
for the Picard solution with bo = (0, 0, —1/2, —1/2) and 

hi{t) =t{t-l)Hi{t) -t/A+l/S (55) 

for the solution with parameters (53). 

Now hi{t) is obtained using the birational canonical transformation (16) with param- 
eters (38) and h{t) replaced with hQ{t). Using (54-55) one can arrive at the following 
answer: 



4dii\z,t) At{t-1) 



TT r, „ , 1 M/ iy sn(z,t)cn(z,t) 



2K{t)' ° ' ^1 /jo; ■ ^^^^ dn{z,t) 

(56) 

In fact, the equation satisfied by ho{t) and hi{t) is the same and it is easy to check that 

71" 71" 

hi{t,x,y) = ho{t,x + -,y + -) (57) 

and 

H,it, X, y) = Ho{t, x+|,?/+|) + - + 4(^3^, (58) 

where we show explicitly a dependence on the fixed parameters x and y. 
Taking (58) into account it is easy to integrate Hi{t) using formula (49). 
The answer is given by the following 

Theorem 4.2 The r-function Ti{t) is given by 

Ti(t) = exp{ J H,{t)dt} = c^ix,y) //-^(l - t)'/' 
where C\ (x, |/) is an arbitrary integration constant. 



(59) 



If we define the sequence of tau- functions T^(i) corresponding to Picard type solutions 
with parameters 

b„ = (0,0,-l/2 + m,-l/2), meZ (60) 
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and two initial conditions 



^ emr) ' '^^"^ ^ ^ ^ ^4(0|r) ' ^^^^ 

then other tau-functions Tjn{t), for m > 1 or m < can be calculated from the difference- 
differential equation 



d 
di 



t(t-l)^logT^(t) 



+ I m - - 



J- m+l 



where c(m) is determined by a normalization of tau-functions. 

Note that the expressions for T„i{t), m ^ 0,1 will be more complicated and involve 
explicitly the function S{u,t) defined in (26). We will not calculate them here. 



5 Conclusion 

In this paper we constructed a sequence of tau-functions for the Picard type solutions of 
Painleve VI equation with parameters (60). In fact, starting with the tau-function for the 
Picard solution (49) one can write many different birational canonical transformations 
and calculate corresponding sequences of tau-functions. We have successfully applied this 
approach to sum up some infinite form factor expansions for the 2D Ising model. All 
details will be given in forthcoming publications. 
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